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Monte Carlo techniques with importance sampling have been extensively applied to 
lattice gauge theory in the Lagrangian formulation. Unfortunately, it is extremely 
difficult to compute the excited states using the conventional Monte Carlo algo- 
rithm. Our recently developed approach: the Monte Carlo Hamiltonian method, 
has been designed to overcome the difficulties of the conventional approach. In this 
paper, we extend the method to many body systems and quantum field theory. The 
Klein-Gordon field theory is used as a testing ground. 

1 Introduction 

Quantum theory is one of the most important achievements in modern science, 
but a lot of non-perturbative aspects are still poorly understood. There are 
two standard formulations in quantum theory: Hamiltonian and Lagrangian. 
A comparison of the conventional approaches is given in Tab.0. 

The Lagrangian formulation is very suitable for applying the Monte Carlo 
(MC) method to systems with many degrees of freedom, and. in the last two 
decades, it has been widely applied to lattice gauge theoryiirBtl In the stan- 
dard Lagrangian MC method, however, it is extremely difficult to compute 
the spectrum and wave function beyond the ground state. On the other hand, 
the standard Hamiltonian formulation is capable of doing it. If fact, some ef- 
forts have been made to solve the lattice QCD Schrodinger equationn. But 
due to the complexity of non-Abelian gauge theory and high dagrees of free- 
dom, only the ground stateO and the lowest lying glueball statesErQtl have been 
investigated. There have also been some attempts to perform numerical sim- 
ulations in the Hamiltonian formulation,. e.g the quantum MC techniqudj, 
MC method using a guided random walkEjEJ, and the MC Green's function 
methodt3ll3. For an review^jf MC methods applied to solve the quantum 
many-body problem see Ref.tfl. 
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Table 1. Comparison of the conventional methods in the standard formulations. 



Formulation 


Hamiltonian 


Lagrangian 


Approach 


Schrodinger Eq. 
H\E n ) = E n \E n ) 


Path Integral 

, , f\d<p]0\w] exp(-SM/fi) 
^ ' fld V ]exp(-S[tp]/K) 


Algorithm 


Series expansion, 

variational, 
Runge-Kutta ... 


MC simulation 


Advantage 


Both the ground state, 
and the excited states 
can be computed. 


It generates the most 

important configs. 
for the measurements. 


Disadvantage 


Analytical methods 
are too tedious for 
many body systems; 
Runge-Kutta works 
only in 1-D. 


It is difficult to study 
the excited states, and 
finite density QCD. 



A natural question is whether in Lagrangian MC simulations one can 
construct an effective Hamiltonian? If yes, the excited .slates can also 
be computed. We have recently proposed a new approacht3 (we called it 
Mnntp rarln HaiTiiltnTiiari method) to investigate this problem. A lot of 
models^'^'EZni3'Ej'E3 in quantum mechanics (QM) have been used to test 
the method, as briefly reviewed in Sect. Q In this paper, we will extend the 
method to many body systems and quantum field theory (QFT). 

2 Algorithm 

2.1 Effective Hamiltonian 

The basic idea of the MC Hamiltonian method for one body quantum mechan- 
ics has been described in details in Rcf.ll3. The (imaginary time) transition 
amplitude between an initial state at position Xi, and time U, and final state 
at Xf and tf is related to the Hamiltonian H by 

M fi = {x f ,t f \x % ,U) = (xfle-^f-^Xi) 

OO 

= J2(xf\Zn)e- E " T ' n (E n \x i ), (1) 

n=l 
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where T = tt — According to Feynman's path integral formulation of QM, 
the transition amplitude is also related to the path integral: 

M fi = J[dx]eMS[x]/h)\ x jJf, (2) 

where S = So + Sy is the action for a given path, So = f t f dt mx 2 /2, and 
Sv=fl f dt V{x). 

The results above-are standard ones in Lagrangian formulations. What 
is new in our methoduJ is to construct an effective Hamiltonian H e s (finite 
N x N matrix) through 

N 

M fi = {x f \e- H ^\ Xi ) = J2(xf\K S )e- E * !T/h (Ef\xi)- (3) 

n=l 

The eigenvalues £^ ff and wave function |££ ff ) can be obtained, by diag- 
onalizing M using a unitary transformation M = WDU, where D = 
diag(e~ El T / h , ...,e~ E]s > T / h ). Once the spectrum and wave functions are 
available, all physical information can also be obtained. 

Since the theory described by H, which has an infinite basis in Hilbert 
space, is now approximated by a theory described by a finite matrix H c g, 
which has a finite basis, the physics of H and H c g might be quite differ- 
ent at high energy. Therefore we expect that we can only reproduce the 
low PT1p rgy jPliysks. of the system. This is good enough for our purpose. 
In RefsMlilijMyy, we investigated many 1-D, 2-D and 3-D QM models 
(Tab.||) using this MC Hamiltonian algorithm. We computed the spectrum, 
wave functions and some thcrmodynamical observables. The results are in 
very good agreement with those from analytical and/or Runge-Kutta meth- 
ods. 



2.2 Basis in Hilbert Space 

To get the correct scale for the spectrum, the position state \x n ) at initial 
time ti or final time tt (Bargman states or box states) should be properly 
normalized. We denote a normalized basis of Hilbert states as |e„), n = 
1, ...,N. In position space, it can be expressed as 

/ \ _ r V >/Ax n , X G [x n , Xn+i] / .s 

0, x<t[x n ,x n+1 ] W 

where Ax n = x n +i — x n . 

There is some arbitrariness in choosing a basis for the initial and final 
states. The simplest choice is a basis with Ax n = const., which is called the 
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Table 2. QM systems, investigated by the MC Hamiltonian method with a regular basis. 



System 


Potential 


QM in 1-D 


V(x) = 
V{x) — ^muj 2 x 2 
V(x) = -F sech 2 (x) 
V{x) = \x 2 + ±x 4 
V(x) = ±\x\ 

oo, x < 
V {X) i Fx 1 x>0 


QM in 2-D 


V(x, y) — ±moj 2 x 2 + ^muj 2 y 2 
V(x, y) — ^muj 2 x 2 + ^moj 2 y 2 + Xxy 


QM in 3-D 


V(x, y, z) = jmL0 2 x 2 + ^muj 2 y' 2 + ±mu} 2 z 2 



"regular basis" . In Refs.E£Hl3't3'E3, the regular basis was used. For many body 
systems or QFT, the regular basis will encounter problems. For example, in a 
system with a 1-D chain of oscillators (see later), if the number of oscillators 
is 30, the minimum non-trivial regular basis is N = 2 30 = 1073741824 w 10 9 , 
which is prohibitively large for numerical calculations of the .matrix elements. 

Guided by the idea of importance sampling, in RefsJ 19 H 2C l, we proposed 
to select a basis from the Boltzmann weight proportional to the transition 
amplitude between x\ — at t[ = and x 't — x n at some t'* . In a free particle 
or harmonic oscillator case, the distribution is just a Gaussian 



where a = y/hT'/m for the free case and a = y/ hsinh(ujT') / (muj) for the 
harmonic oscillator. We call such a basis the "stochastic basis", to be used 
as states at initial time U or final time tf. For more complicated models, one 
may still use the Gaussian distribution, but with an adjustable parameter a. 

2.3 Matrix elements 

As explained above, the calculation of the transition matrix elements is an 
essential ingredient of our method. The matrix element in the normalized 
basis is related to {x n i, tf\x n , ti) by 

M n , n = (e n ,,t f \e n ,ti) = r >+1 dx 1 f n+1 dx" < x '^f\ x "'^ 
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Paths of a free particle generated by MC 
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Figure 1. Config. stands for the classical path, Config. 1, Config. 2 and Config. 3 stand 
the paths generated by MC. Here Nt = 20 (there are 19 time slices between t; and tj). 



yj Ax„iAx n (x n > ,tf\x n ,ti 



(6) 



where for the regular basis, Ax n — const., and for a stochastic basis, Ax n — 
1/ (Pbasis[%n]N) . (x n i , t f\x n , U) can be calculated by MC as follows. 

(a) Discretize the continuous time with time interval At = T/N t . 

(b) Generate free configurations [x] between t £ (U,tf) obeying the Boltzmann 
distribution 



exp(-S [x]/?t) 
J[dx] exp(— So[x]/h) 



\x n ,U ■ 



(7) 



Figure [l] shows a sample of free configurations generated by MC method with 
importance sampling, 
(c) Measure 



(O v ) = I [dx] exp(-Sv(x)/h)\Zl'£' P [x]. 



(8) 
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The path integral in Eq. (0) is then 



(xn',tf\x„,ti) = {O v )\l ^-^exp --^^(x nl - x n ) . (9) 



2hT' 



3 Quantum Field Theory 

3.1 The Klein- Gordon model 

The main purpose of the algorithm is to study many body systems and QFT. 
As an example, we consider a free scalar field theory: the Klein-Gordon model 
in 1+1 dimensions, which has the Hamiltonian 

h = j dx \L^ + ^(^j + n§^. (io) 

Discretizing the theory on a finite string with N osc sites and lattice spacing 
a, and making the replacement x — > Xj — ja, n(x) — > a^ 1 ' 2 H(xj), <p(x) — > 
a _1 / 2 </)(xj), J — ► QsX^fci" an£ i v ~ * a ^-> ^ ne Hamiltonian becomes 

JV„sc 



H = T,\ ( n2 ( x ^ + " 2 (<^) - <K^-+i)) 2 + & 4> 2 (xj)) , (ii) 

which describes a chain of N osc coupled oscillators0 with m = 1 in 1 spatial 
dimension. The field 4>(xj) and its conjugate momentum II(xj) = 4>(xj) at 
the j-th site correspond respectively to the displacement and velocity of the 
j-th oscillator. 

In QFT, the initial and final states are defined by 

= \(f)(xi,ti),...,<t>(x Naac ,ti)), 

\*f,tf) = \<Kx u t f ),...,<l>(x NM ,tf)). (12) 
The transition amplitude between them is 

M fi = ($/,*/!$<, ti) = (^le-^*/-**)/^*,) 

= fme*p(-swh)\%;t' 

,N t -l N OB0 

W«->oo. 



k=l j=l 



i^/ n ff v2^Ai #fe ' tfc) (i3) 
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where the Euclidean action is 

A N t -lN osc 

k=i i=i 



((j>(xj,t k+1 ) - (j>(xj,t k )y 

(AO 2 



+ SI 2 (<t>(xj, tk) - <f>(xj+i, tk)) +VLl4> 2 {x J ,tk) 



(14) 



For a finite -/V osc and non-zero f2, one has to introduce a boundary condition 
for ^>(:x;jv 03C +i)- We implement the periodicJioundary condition 4>{ x N oac +x) = 
4>(xi). The spectrum is analytically knowncd: 



E 



(15) 



where oji = y fig + 4£1 2 sin (piAx/2), and ni, ...,njv osc =0,1,... . We can 
also compute the transition amplitude analytically using the discrete Fourier 
transformation 

Ap 



<t>{ x j>t) = y2-ZF=wp(ipiXj)4>(pi,t), 
, V 2ir 



Ax 



V Z7T 



(16) 



where Ax = a, Ap = 2n/ (N osc Ax), x 3 = [-{N osc - l)/2 + (j - l)]Ax, 
pi = [-(N osc - l)/2 + (I - l)]Ap, and j, I = 1, .., iV osc . The result is 

JVo,c 



M fi = n 



i=i 



moji 



2Trhsmh(uJiT) 



exp 



ojiAp/ Ax 
2?ismh(ujiT) 



x {(\kputf)\ 2 + \kpuU)?) cosh(«,r) - 2Re ^(p,, ^^(pj.f,) 



(17) 



Stochastic basis 

In QFT, a stochastic basis has to be used for the initial and final states, 
because in the thermodynamical limit, N osc — > oo. Although the variables 
are coupled in position space, we may still use a stochastic basis [$ n ], 
n = 1, N according to the distribution with an adjustable parameter a 

N 



Pi 



bo 



sis Y^n 



n4 



27TCT 



exp 



</£(sj) 
2a 2 



(18) 
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The matrix element between the normalized initial and final states is 



tf\e n ,ti) « y/A$ n ,A$ n {$ n ,,tf\$ n ,ti), (19) 
tesimal volumes for the initic 

A$„ = d(f> n (xi)...d<f> n (x Noso ) 



where the infinitesimal volumes for the initial and final states are 

1 



Pbas^sl<5>n}N , 

A$„/ = d<p n ,{x 1 )...d^ n ,{x Nosc ) = fj— 4jO]v - ( 2 °) 



4 Results 



We consider N osc = 9, a = 1, O = 1, Oq = 2, m = 1, H = 1, and T = 
2. For the adjustable parameter c in the stochastic basis, we choose a = 
hsmh(QoT') / (mQo) with T" = T for simplicity. After the stochastic basis 
with ./V = 1000 is generated, we obtain the matrix elements M n > n using Eqs. 
(|l7|), (|Ts| ) and (|20|). (In the future work, we will calculate (<& n i ,tf\$ n ,ti) 
directly by MC). Then we co mpu te the eigenvalues and eigenvectors using 
the method described in Sect. [O . 

Table | gives a comparison between the spectrum from the effective 
Hamiltonian with the stochastic basis and the analytical formular Eq. ( |l5| ) 
for the first 20 states. They agree very well. This means that an arbitrary 
choice of Gaussian distribution for the stochastic basis is good enough. (Of 
course, one should study systematically the dependence of the results on a) . 

We have also computed thermodynamical quantities such as the partition 
function average energy E and specific heat C . The analytical results are 

Z{(3) = Tr(exp(-/?£0) = [] 



i= 



L 2sinh(/3nw;/2) : 



£708) = |Tr (ffcxp (-(3H)) = -^5§^ = E IT coth (^/ 2 ) . 

^.^..^.^K-^-)-, (21) 

where (3 = T/h, the temperature T = l/(/3fcs), and &b is the Boltzmann 
constant. Since we have approximated H by H c g , we can express those ther- 
modynamical observables via the eigenvalues of the effective Hamiltonian 

JV 



luo 'mch: submitted to World Scientific on February 1, 2008 8 



Table 3. Comparison of the spectrum of the Klein-Gordon model on the lattice, between 
the MC Hamiltonian method with a stochastic basis and the analytic ones. 



71 




^cxact 


1 


10 


904fifi31921fi8 


10 


944060480668 


2 


12 


956830557334 


12 


944060480668 


3 


12 


985023578737 


13 


057803869484 


4 


1 3 


04431 1 582647 


1 3 


057803869484 


5 


1 3 


299967341 242 


1 3 


321 601 993380 


6 


13 


345480638394 


13 


321601993380 


7 


13 


552195133687 


13 


589811791733 


8 


13 


585794986361 


13 


589811791733 


9 


13 


680136748933 


13 


751084748745 


10 


13 


744919087477 


13 


751084748745 


11 


14 


984737011385 


14 


944060480668 


12 


15 


012353803145 


15 


057803869484 


13 


15 


057295761044 


15 


057803869484 


14 


15 


108904652020 


15 


171547258300 


15 


15 


125356713561 


15 


171547258300 


16 


15 


187413290039 


15 


171547258300 


17 


15 


308536490102 


15 


321601993380 


18 


15 


396255686587 


15 


321601993380 


19 


15 


420708031412 


15 


435345382196 


20 


15 


432823810789 


15 


435345382196 



N 



c-w = ^ (x; (£ ^y - (^csjj'j • (22) 

Since this is a static system, the eigenvalues should in principle not vary with 
(3, which is also confirmed numerically within statistical errors when (3 and 
N are not too small. Using this assumption and the spectrum at j3 = 2, 
we obtain the thermodynamical quantities for other (3. The results for the 
free energy, average energy and specific heat as a function of (3 are shown in 
Figs. ^, ^ and |]. Again, the results from the MC Hamiltonian are in good 
agreement with the analytical ones when [3 > 1. 

Preliminary results for larger N osc indicate that it is not necessary to 
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Figure 2. Free energy of the Klein-Gordon model on a 1+1 dimensional lattice. 



increase N accordingly. This property is very important for a feasible appli- 
cation of the algorithm to many body systems and QFT. 

5 Summary 

In this paper, we have extended the effective Hamiltonian method with a 
stochastic basis to QFT, and taken the Klein-Gordon model as an exam- 
ple. The results are very encouraging. We believe that the application of 
the algorithm to more complicated systems will be very interesting for non- 
perturbative investigation beyond the ground state. 
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Figure 3. Average energy of the Klein-Gordon model on a 1+1 dimensional lattice. 
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